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The coupling of distinct systems underlies nearly all physical phenomena. A basic instance is 
that of interacting harmonic oscillators, giving rise to, for example, the phonon eigenmodes in 
a lattice. Of particular importance are the interactions in hybrid quantum systems, which can 
combine the benefits of each part in quantum technologies. Here we investigate a hybrid 
optomechanical system having three degrees of freedom, consisting of a microwave cavity 
and two micromechanical beams with closely spaced frequencies around 32 MHz and no 
direct interaction. We record the first evidence of tripartite optomechanical mixing, implying 
that the eigenmodes are combinations of one photonic and two phononic modes. We identify 
an asymmetric dark mode having a long lifetime. Simultaneously, we operate the nearly 
macroscopic mechanical modes close to the motional quantum ground state, down to 1.8 
thermal quanta, achieved by back-action cooling. These results constitute an important 
advance towards engineering of entangled motional states. 
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One of the present goals in physics is the explanation of 
macroscopic phenomena as emerging from the quantum- 
mechanical laws governing nature on a microscopic scale. 
This understanding is also important for future quantum 
information applications^'^, as many of the most promising 
platforms are based on nearly macroscopic degrees of freedom. 
For macroscopic mechanical objects, their potential quantum 
behaviour^ has been actively investigated with resonators 
interacting with an electromagnetic cavity mode^. In particular, 
freezing of the mechanical Brownian motion^'^ down to below 
one quantum of energy has recently been observed^'^^. 

Another result of radiation-pressure interaction is the mixing 
of the normal modes of vibration into linear combinations of the 
uncoupled phonon and photon modes. This has been verified 
with optomechanical Fabry-Perot cavity^ and by its on-chip 
microwave analogue^ ^ using an aluminium membrane. All this 
work is paving the way towards engineering non-classical 
motional^ and hvbrid quantum states ^'^'^'^^ for basic tests of 
quantum theory as well as applications in foreseeable 
future. 

Corresponding phenomena in systems comprising more than 
two active degrees of freedom^^"^ have received less attention. 
In the optomechanical crystal set-up, coupling of mechanical 
vibrations via radiation-pressure interaction has been demon- 
strated with the zipper cavity^^, however, with some direct 
interaction between the beams. In the microwave regime, such 
measurements have remained outstanding. In this work, we take 
a step further, and examine a multimode system where 
two micromechanical beams, having resonant frequencies 
«i/(27i) = 32.1 MHz and cl)2/(27i) = 32.5 MHz, are each 
coupled to a microwave on-chip cavity. We obtain the first 
evidence of hybridization of all the three degrees of freedom. 
This was made possible by operating, as opposed to the optical 
set-up, in the good-cavity limit, that is, both coi and CO2 were 
much larger than the cavity decay rate 7c- Simultaneously, the 
mechanical modes are found at low occupation numbers n^y 
near the quantum ground state, down to = 1.8 ± 0.5. These 
are the lowest occupations recorded with nanowire resonators to 
date, and they occur via the sideband cooling, starting from 
operation at dilution refrigerator temperatures at a few tens of 
millikelvin. 



Results 

Cavity with several mechanical modes. In the optical domain^^, 
a three-mode optomechanical system can be pictured as a Fabry- 
Perot cavity with both massive mirrors mounted via springs, see 
Fig. la. In the microwave version, (Fig. lb), the mirrors are 
replaced by movable capacitors formed by the mechanical reso- 
nators. The electromechanical coupling arises when each 
mechanical resonator (labelled with j), expressed as time- varying 
capacitors Cp independently modulates the total capacitance, and 
hence, the cavity frequency cOc- This is described by the coupling 
energies gj — {(Dj2C)dCjldxj. In the actual device (see the 
micrograph in Fig. Ic), the mechanical resonators are spatially 
separated by about 100 jim, and we hence expect the direct 
interaction between the beams via the substrate material to be 
negligible. We drive the cavity strongly by a microwave pump 
signal applied at the frequency cop near the cavity frequency. The 
pump induces a large cavity field with the number of pump 
photons np ^ 1. This allows for a linearized description of the 
electromechanical interaction and results in substantially 
enhanced effective couplings Gj—gjXQjy/n^, where 
Xqj — ^h/2mjC0j is the zero-point amplitude. 

The general Hamiltonian for N mechanical resonators with 
frequencies (Op individually coupled to a common cavity mode. 
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Figure 1 | Cavity-controlled coupling of micromechanical resonators. 

(a) Optomechanical system with both mirrors movable has three degrees 
of freedom. These are the optical field inside the Fabry-Perot cavity and the 
displacements x-\ and X2 of the mirrors, (b) Analogue of panel a with a 
microwave-regime cavity represented by the lumped element inductors and 
capacitors and arbitrary number y = 1 ... N of mechanical resonators. The 
system is probed through the coupling capacitor Q by shining a strong 
pump microwave tone (frequency cop), which corresponds to the laser in a, 
and possibly a weak probe tone (frequency cojn and voltage amplitude Vjn) 
in a dilution refrigerator down to 25 mK temperature. The outgoing wave 
Vout can be monitored, (c) Images of the all-aluminium superconducting 
sample, showing four fabrication jigs for a beam. Total of three beams (x-\ ... 
X3) were operating. Beams 1 and 2 had large electromechanical couplings 
gi/27r = 1.8MHznm-\ ^2/271 = 2.0 MHz nm whereas the third beam 
had an order of magnitude smaller coupling. The scale bars are 50|im on 
the left and 10|im on the right. 



written in a frame rotating with the pump, becomes 

H= -HAa^a^H'^Wjbpj-H{a'^ ^a)'^Gj\^b^ ^bj] (1) 

j=i j=i 

Here, a is the annihilation operator for the cavity mode, while bj 
are those of the mechanical resonators. If one assumes roughly 
similar mechanical resonators, viz. coj ~ co^ and Gj ^ G, the 
cavity becomes nearly resonant to all of them at the red- sideband 
detuned pump condition A = — co^. 

Coupling of the mechanical resonators via the cavity 'bus' can 
be anticipated to be significant if the mechanical spectra begin to 
overlap when their width increased by the radiation pressure, 
7eff = 7m + 7opt> grows comparable to the mechanical frequency 
spacing. Here, 7opt = 4G^/7c- The equations of motion following 
from equation (1) allow one to vehfy the above assumption. We 
will in the following focus on two resonators, N=2. The result of 
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such a calculation, which is detailed in the Supplementary 
Methods, is that they experience an effective coupling with energy 
Gi2 = GiG2/7c- 

In the limit of strong coupling, ^/(Cl+~G^y2 ;:$> y^/A, the 
system is best described as a combination of two 'bright' modes 
with linewidths ~ 7^-/2, and of a 'dark' mode having a long 
lifetime ~ 7^. In the latter, the two mechanical resonators 
oscillate out-of-phase and synchronize into a common frequency 
(«i + a92)/2. In any case, the mode structure can be written with 
the normal coordinates for mode n, for position quadrature 
according to X*^"^ ^Dc^^'^c + X^jLi^j"^^;- Here we used the 
dimensionless amplitude of the mechanical oscillations 

= (I/V2) (^J + ^;), and similarly the dimensionless flux 
qc=[l/V^ W+(i) in the cavity. The mode expansion coeffi- 
cients give the participation ratio of each uncoupled mode. They 
satisfy (dW)^+E7.,(<)^ = 1. 

In the experiment, we use doubly clamped flexural beams, 
which are made by the use of ion beam milling of aluminium^ \ 
having an ultranarrow ^ lOnm vacuum slit to the opposite 
end of the cavity for maximizing the coupling. The cavity is a 
half-wavelength floating microstrip similar to the one used 
in Massel et al}^, resonating at CDj{2n) = 6.98 GHz. The total 
cavity linewidth 7c = 7e + 7i — (27i) x6.2 MHz is a sum of the 
internal damping 7i/(27i) = 1.4MHz, and the external damping 
7e/(27i) = 4.8 MHz. At the highest powers discussed, however, we 
obtain a decreased 7i/(27i) = 0.9MHz, typical of dielectric loss 
mechanism. 

There are a total of three beams as shown in Fig. Ic. Two of the 
beams have a large zero-point coupling of giXoi/2n — 39 Hz and 
^2^02/271 = 44 Hz. The frequencies of beams 1 and 2 are relatively 
close to each other, cl)2 — ft^i (27i)450kHz, such that it is 
straightforward to obtain an effective coupling G12 of the same 
order. The third beam has an order of magnitude smaller 
coupling, and hence we neglect it in the full calculation, setting 
N=2 in equation (1) in what follows. 

The measurements are conducted in a cryogenic temperature 
below the superconducting transition temperature of the Al film 



1 K) to reduce losses. We used a dilution refrigerator set-up as 
in Massel et al}^ down to 25 mK, which allows for a low 
mechanical mode temperature in equilibrium. Apart from this, 
cryogenic temperatures or superconductivity are not necessary^^. 

Experimental data. The modes can be probed by adding to the 
experimental set-up another, the probe, tone^^'^^'^ with 
frequency cOin- The maximum pump power we can reach is 
4 jiW, inducing a coherent np ^ 6x10^ photon occupancy in 
the cavity, and effective coupling G12 = (27i) x 150 kHz (with 
Gi = 0.93 MHz, G2= l.OMHz). In Fig. 2a,b, we display changes 
of the cavity absorption by increasing the pump power. The two 
peaks at the bottom of the cavity dip are due to beams 1 and 2. 
They broaden with increasing np, finally leaving behind a narrow 
dip between them. The overall lineshape is clearly not a sum of 
two Lorentzian peaks. 

In Fig. 2c, we show a zoom-in view of the two peaks, together 
with theoretical predictions. One immediately sees that an attempt 
to simulate either peak by a single mechanical resonator coupled to 
cavity {N—\ in the analysis), thereby neglecting the cavity- medi- 
ated coupling, fails to explain the lineshape at large np > 10^. The 
full simulation with both beams 1 and 2 included, however, pro- 
duces a remarkable agreement to the experiment. To create these 
curves, we further take into account a shift of the mechanical fre- 
quencies due to second-order interaction^, given as — {\/2)g-^^ n^, 

where gf^ = xI.(Dc/2C {d^q/dxf^ - (27r)3 x 10 " ^ Hz is the 

second-order coupling coefficient. The third beam, visible in 
Fig. 2c as a narrow feature just below 6.982 GHz, was later fitted by 
running a single-resonator calculation. This is justified because its 
mixing is negligible owing to weak coupling. 

The narrow dip between the peaks manifests the onset of the 
dark mode, where the cavity participation approaches zero with 
growing Gp as occurs also in other tripartite systems^^. As usual for 
weakly decaying states, the dark modes can be usefijl for storage of 
quantum information^^'^^. The other two dips (best discerned in 
the theoretical plots at higher Gp see Fig. 2d) are the bright modes. 
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Figure 2 | Hybridization of microwave photons and two radio-frequency phonons. (a) Measured reflection coefficient Sn = Vout/^n of the probe 
microwave while a strong pump is applied at the red sideband, A ~ —coi, at increasing pump amplitudes, from bottom to top. A low value of S^^ means 
high resonance absorption. (b,c) Zoomed-in view of panel a. An excellent agreement with the theory is obtained with reflection calculation using a 
complete model based on equation (1), for two beams 1 and 2 (solid lines). The dashed lines show a fit with an independent-resonator model for beam 1. 
The narrow feature at 6.982 GHz is due to beam 3. The curves are displaced vertically by 3dB for clarity, (d) Theory plot for higher np; the dark mode is 
expected to become narrower and more and more phonon-like. Here, A= — {(0^ + 0)2) /2. The other parameter values are G/yc = 0.28 and 0.63, 
np = 2 X 10^ and 10^^, for the lower and upper curve, respectively, (e) Calculated dark mode decomposed into a linear combination of a photon and two 
phonons. Blue denotes the cavity component D^, while red (D^) and black (D2) are components of the mechanical resonators 1 and 2, respectively. The 
parameters are as in panel c. The arrow denotes the highest experimental value. For simplicity, we take gj^^ =0 in panels d and e. 
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and they retain a fully mixed character of all three degrees of 
freedom. In Fig. 2e, we display the prediction for the mode 
expansion coefficients with respect to coupling. Asymmetry in the 
role of the mechanical modes is sensitive to parameters. With the 
maximum coupling in the experiment (Gj/y^ ~ 0.15), the normal 
modes are well mixed: the dark mode can be expressed as 
[De,Di,D2] - [0.3,0.3,0.9]. 

One can also vary the pump detuning, thereby effectively 
detuning the cavity and the mechanical resonators^ An 
anticrossing of the cavity and mechanical modes is seen in 
Fig. 3a up to detuning A — cOm, above which the system exhibits 
parametric oscillations. The pertinent simulation. Fig. 3b, portrays 




6.965 
6.935 



Figure 3 | Pump-detuning measurements, (a) Measured spectrum of the 
tripartite cavity electromechanical system when pumped at a fixed input 
power such that at resonance, A ~ ( — coi, — 0J2), — 5.6x10^. The pump 
frequency was slowly swept from left to right, (b) Theory prediction. 



the main features of the measurement, except bending of the cavity 
frequency towards lower value in an upward sweep of A. We expect 
such possibly hysteretic behaviour of the cavity frequency to be due 
to second-order effects beyond the present linear model. 



Sideband cooling. We now turn the discussion into showing how 
the tripartite system resides nearly in a pure quantum state during 
the mode-mixing experiments. The mechanical resonators can be 
sideband cooled with the radiation pressure back-action^'^^"^^ 
under the effect of the pump, which also creates the mode mixing. 
Although the ground state of one mode was recently achieved this 
way^'^^, cooling in multimode or nanowire systems is less 
advanced. The efficiency of cooling with two mechanical 
resonators is supposed to become hindered at those values of 
effective coupling where mode mixing starts to set in^^. At modest 
pump powers, relevant for achieving the lowest phonon numbers 
in the present case, mixing remains insignificant. 

The incoherent output spectrum can be due to either a finite 
phonon number n^^ or non- equilibrium population number nj 
of the cavity. The former manifests as a Lorentzian centred at the 
upper sideband, as in Fig. 4a, with the linewidth 7eff = 7m + 4G^/ 
7c, whereas the latter gives a small emission at broader bandwidth. 
The analytical approximation^ and experimental details are 
discussed in Methods. 

At low input power, when the cavity back-action damping is 
small and nj ~ 0, phonon number is set by thermal equipartition 



n^hcDj — k^T. This is best observed by the area of the mechanical 
peak in the output spectrum, as in Fig. 4a, which should be linear 
in cryostat temperature. The linearity holds down to about 
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Figure 4 | Sideband cooling of motion down to 1.8 thermal quanta, (a) Cavity output spectra representing the thermal motion peaks of beams 1 (left) and 
2 (right) measured at 345 mK (red) and 35 mK (black), with np = 4x10^. (b) Calibration of the detection according to temperature dependence of the 
mechanical energy rimHcom^kBT, given by the area of the peak (data for beam 2). The straight line is a linear fit to the data above 120 mK. (c) Output 
spectra of beam 2 at pump occupation np, increasing from top to bottom. The curves are shifted vertically by 0.3 y-axis units for clarity. The corresponding 
mechanical occupations are written on the right, (d) Progress of back-action cooling (beam 2) as a function of pump occupation. The mechanical mode 
temperature is given by circles, and the cavity temperature by squares. Red: high temperature, 245 mK; blue, base temperature (28 mK). The error bars are 
from the fits in panel c, given as one s.d. The ideal theoretical behaviour are plotted by the solid lines. Unless the data point for nj, including its error bars, 
was of positive sign, only the absolute value of the error bar is plotted. 
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150 mK temperatures (Fig. 4b), below which we observe 
intermittent heating, which is sensitive to parameters. 

The cooling of beam 2 is observed in the mechanical peaks 
shown in Fig. 4c (see Methods), and the resulting phonon 
numbers are displayed in Fig. 4d. At a relatively high temperature 
of 245 mK, the data points follow theory well up to input powers 
corresponding to np'^5xl0^. At the minimum cryostat 
temperature of 28 mK, the mode is not thermalized in a 
manner that depends irregularly on power. We obtain a 
minimum phonon number in beam 2 of nm ~ 1.8, where the 
mechanical mode spends one-third of the time in the quantum 
ground state. The other mechanical mode, beam 1, cooled 
simultaneously down to 2.5 quanta, possibly compromised by the 
slightly smaller coupling. The cooling is quite clearly bottlenecked 
by heating, by the pump microwave, of the bath to which the 
mechanical mode is coupled. For the optimum data points, the 
starting temperature for cooling was raised to 150 mK, and above 
this it grows roughly quadratically with rip. 



Discussion 

The coupling of micromechanical resonators mediated via 
microwave photons in an on-chip cavity is a basic demonstration 
of the control of a multimode mechanical system near the 
quantum limit. The set-up provides a flexible platform for 
creation and studies of non- classical motional states entangled 
over the chip^^, or over macroscopic distances^^'^^. The set-up 
is easily extended to embody nearly arbitrarily many mechanical 
resonators, hence allowing for designing an electromechanical 
metamaterial with microwave- tunable properties. For future 
quantum technology applications even at elevated tempera- 
tures^^, data may be transferred for storage^^ from photonic 
modes into the dark states, which decay very slowly at the 
mechanical damping rate. 



Metliods 

Device fabrication. The lithography to define the meandering cavity and jigs for 
the beams consisted of a single layer of electron-beam exposure, followed by 
evaporation of 150 nm of aluminium on a fused silica substrate. As the cavity is of 
planar structure, it has no potentially problematic cross-overs, which tend to have 
weak spots in superconducting properties and hence could heat up the cavity. We 
thus expect the cavity to portray a high critical current limited only by the intrinsic 
properties of the strip. With a 2 |im wide strip, we expect the critical current to be 
several milliamperes. With the maximum circulating currents in cavity (Fig. 2), 
about 0.5 mA of the peak value, there was no sign of non-linearity in the cavity 
response. 

The mechanical resonators were defined by focussed ion beam cutting, as in 
Sulkko et al?^. We used low gallium ion currents of 1.5 pA, which gives the 
nominal beam width of about 7 nm. To maximize sputtering yield with minimal 
gallium contamination, we used a single-cutting pass mode. Otherwise, the cut 
beams tend to collapse to the gate. With the current recipe, fabrication of down to 
10 nm gap widths over 5-10 |am distance can be done with about 50% yield. 

Theory of cavity-coupled resonators. The Hamiltonian for N mechanical reso- 
nators each coupled to one cavity mode via the radiation-pressure interaction is 

;=1 ;=1 



The pump with the Hamiltonian Hp = iH^P^y^/HcOcia'^ + a) cos co^t induces a 
large cavity field with the number of pump photons ^ I. 

Under the strong driving, the cavity-mechanics interaction can be linearized 
individually for each beam, resulting in equation (1). Further details of the cal- 
culations are given in the Supplementary Methods. 

Sideband cooling. For one beam, the measured output spectrum divided by sys- 
tem gain, Sout(<^)> carries information on the phonon number according to^ 



rim - 2n, 



(3) 



which is a Lorentzian centred at the cavity frequency. A non-zero base level is due 
to possible broadband emission from the cavity, due to a thermal state with 



occupation nj. At modest pump (cooling) powers of utmost interest, coupling of 
the mechanical resonators remains insignificant, and thus we can apply 
equation (3) separately for each resonator. Some deviations can, however, occur at 
the highest powers with ~ 3x10^. 

To obtain the phonon number under strong back-action, we fit a Lorentzian to 
each peak as in Fig. 4c, obtaining independently the base level, amplitude and 
linewidth for every input power. These values are then compared with equation (3). 
This leaves, however, yet too many unknowns to obtain the occupations of the 
mechanics and cavity. This is basically because of uncertainties in the attenuations 
of both input and amplifier lines, which limit the accuracy of estimating rip. We can 
obtain a calibration using the linear temperature dependence, and by fitting the 
dependence of yeff on input power. 

As opposed to, for instance, Rocheleau et al.^, the bottleneck for cooling is not 
heating of the cavity by strong pump, as seen in Fig. 4d, where the cavity tem- 
perature is increasing only little up to the optimum cooling powers corresponding 
to about 0.15 mA of peak current values in cavity. This can be due to the simplistic 
structure of the cavity. 

The above analysis supposes the presence of nearly no non-equilibrium photons 
in the cavity at zero or low input powers. This situation can be investigated by 
observing a possible emission from about cavity linewidth under these conditions. 
A complication arises because the cavity linewidth is so large that a small change is 
easily overwhelmed by modest standing waves in cabling. We avoided this issue by 
using temperature dependence of the cavity frequency, due to the kinetic induc- 
tance in the long microstrip. Thereby, a broad emission peak moving according to 
the cavity frequency would be easily distinguishable. We observed no such signal 
down to the level of nj<0.1, which justifies the above analysis. 
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